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Abstract — This paper shows how to reduce the otherwise hard 
joint relay positioning and flow optimization problem into a 
sequence a two simpler decoupled problems. We consider a 
class of wireless multicast hypergraphs mainly characterized by 
their hyperarc rate functions, that are increasing and convex in 
power, and decreasing in distance between the transmit node 
and the farthest end node of the hyperarc. The set-up consists 
of a single multicast flow session involving a source, multiple 
destinations and a relay that can be positioned freely. The 
first problem formulates the relay positioning problem in a 
purely geometric sense, and once the optimal relay position is 
obtained the second problem addresses the How optimization. 
Furthermore, we present simple and efficient algorithms to solve 
these problems. 

I. INTRODUCTION 

We consider a version of network planning problem under 
a relatively simple construct of a single session consisting of 
a source s, a destination set T and an arbitrarily positionable 
relay r, all on a 2-D Euclidean plane. The problem can then 
be stated as: What is optimal relay position that maximizes 
the multicast flow from s to T? Similarly, we can also ask: 
What is the optimal relay position that minimizes the cost (in 
terms of total network power) for a target multicast flow F? 

A fairly general class of acyclic hypergraphs are considered. 
The hypergraph model is characterized by the following rules 
of construction of the hypergraph Q{I\f,A): 

1) GiM, A) consists of finite set of nodes TV positioned on 
on a 2-D Euclidean plane and a finite set of hyperarcs A. 

2) Each hyperarc in A emanates from a transmit node and 
connects a set of receivers (or end nodes) in the system. 
Also, each hyperarc is associated with a rate function 
that is convex and increasing in transmit node power and 
decreasing in distance between the transmit node and the 
farthest node spanned by the hyperarc in the system. 

3) Each end node spanned by the hyperarc can decode the 
information sent over the hyperarc equally reliably, i.e. 
all the end nodes of an hyperarc get equal rate. 

In relation to the special case of our hypergraph model, the 
authors addressed the first question (max-flow) in the context 
of Low-SNR Broadcast Relay Channel in H]. 

This paper has two major contributions. Firstly, we solve 
the general joint relay positioning and max-flow optimization 
problem for our hypergraph model. Secondly, we address the 



min-cost flow problem and establish a relation of duality 
between the max-flow and min-cost problems. An efficient 
algorithm that solves the joint relay positioning and max-flow 
problem is presented, in addition to an algorithm that solves 
an important special case of the min-cost problem. 

The relay positioning problem has been studied in various 
settings Q-l?). In most cases, the problem is either heuris- 
tically solved due to inherent complexity, or approximately 
solved using simpler methods but compromising accuracy. 
We reduce the non-convex joint problem into easily solvable 
sequence of two decoupled problems. The first problem solves 
for optimal relay position in a purely geometric sense with no 
flow optimization involved. Upon obtaining the optimal relay 
position, the second problem addresses the flow optimization. 
The decoupling of the joint problem comes as a consequence 
of the convexity (in power) of hyperarc rate functions. 

The next section develops the wireless network model. Sec- 
tion |lll] presents the key multicast flow concentration ideas for 
max-flow and min-cost flow that are central to the reducibility 
of the joint problem. In Section |IV] we present the algorithms 
and Section |V] contains an example where the results of this 
paper are applied. Finally, we conclude in Section [Vll 

II. PRELIMINARIES AND MODEL 

Consider a wireless network hypergraph Q{MtA) consist- 
ing of \J\f \ = n+2 nodes placed on a 2-D Euclidean plane with 
1^1 number of hyperarcs and the only arbitrarily positionable 
node as the relay r. The node set J\f — {s, r, ti, .., t„} consists 
of a source node s, a relay r and an ordered destination set 
T = {ti, .., t„} (in increasing distance from s). Their positions 
on the 2-D Euclidean plane are denoted by the set of two-tuple 
vector Z = {zi ~ {xj,yj)\ij e N]. 

All hyperarcs in A are denoted by (u, T4,^ ), where u is 
the transmit node and T4„ = {vi,..,Vk^} is the ordered 
set (in increasing distance u) of end nodes of the hyperarc, 
and Vk^ C M\{u]. The hyperarcs emanating from a trans- 
mitter node are constructed in order of increasing distances 
of the receivers from the transmitter (refer Figure [T]i. This 
construction rule captures the distance based approach and is 
analogous to time sharing for broadcasting. Note that, this is 
one technique to construct the hypergraph (A/^, A), our model 
allows arbitrary styles of hypergraph construction that follow 




Fp . Define F to be the the multicast flow 



Fig. 1. Hyperarcs are constructed in increasing order of distance from the 
transmitter, (a)-(c): 3 node system, (a): Source liyperarc set - {{s, r), {s, rt)}. 
(b): Relay liyperarc set - {(r, t)}. (c): Hypergraph g{Af,A). (d)-(f): 4 
node system with T = {ii,i2} such that Dsr < Dsti < Det^ and 
Drti < Drt2- (d) Source hyperarc set - {{s,r), {s,rti), {s,rtit2)}. (e) 
Relay hyperarc set - {{r,ti), {r,tit2)}. (f): Hypergraph Q{Af,A). 

the above three mentioned rules. Although, since time sharing 
is optimal for broadcasting we will stick to this technique as 
the main example in this paper All the nodes in the set Vk^^ 
receive the information transmitted over the hyperarc (w, Vk^) 
equally reliably. Any hyperarc {u, T4„) S ^ is associated with 
a rate function R"^,^^ = /(-P„" ^ , D^v^^ ), where PJ^^ and D^.^^^^ 
denotes the fraction of the total transmit node power allocated 
for the hyperarc and the Euclidean distance between transmit 
node u and the farthest end node Vk^, respectively. 

The hyperarc rate function i?"^ is increasing and convex 
in power P"^ and decreasing in D^vk^ ■ Furthermore, without 
loss of generality, we write the hyperarc rate function into two 
separable functions of power and distance 



9{Pl)-h{Duv,), (1) 



where g : R+ — > R+ is increasing and convex and h : 
R+ — > R+ is increasing. Mainly, we will be concerned with 
the first equation in ([T]i. Moreover, to comply with standard 
physical wireless channel models we assume that 



< 



dh{D.u 



dD,, 



(2) 



V(i^„"^^^ = i^™,J e doxn{P'^^^^,Duv,J- If the functions g 
and h are not differentiable entirely in dom(P"^ 
then Inequality [2] can be rewritten with partial sub-derivatives, 
implying that differentiability is not imperative. 

Denote the convex hull of the nodes in A/^\{r} by C. For a 
given relay position G C, let Li ~ .., } be the set of 
paths from s to a destination ti G T and let L ~ 
be the set of paths from s that span all the destination set T, 
therefore L C Uie[i «] ^i- Moreover, any path in the system 
consists of either a single hyperarc or at most two hyperarcs 
as there are only two transmitters in the system. Let ji and i/ 
denote the total given power of source and relay, respectively, 
and 7 = ^ denote their ratio, where 7 G (0, 00). Denote with 
Fii and Fi the flow over the path P (for j G [l,Ti]) and the 
total flow to the destination U G T, respectively, such that 



Fi - Eje[i,r,] 

from s to the destination set T as the minimum among the 
total flows to each destination, then for a given relay position 
G C the multicast max-flow problem can be written as. 



Maximize 



F = min Fi 

iG[l,ri] 



subject to; P; < ^ Fp , Vi G [1, 



(A) 



(3) 



0<Fpe€{P,D), VjG [l,T,],ViG (4) 

The hyperarc rate constraints and node sum-power constraints 
are denoted by the set €{P, D) in Program (A) for simplicity. 
Program (A) in general is non-convex, as the path flow 
function Fp can be non-convex, e.g. let the path G Li 

be /] = {{s,Vu^,{r,Vk^)}, = {(s, rti), (r, iiis)} in 
Figure W)), then F,. = min(i?^^^ , Rl^^ ). 

Now we define the notion of cost for a given hyperarc rate 
k... ^ f^ rp,. . — , ni/ given by the 



^Vk^ ~ h{D. 



> 0. The cost of rate i?" 



total power consumed by the hyperarc to achieve K^^ 



(5) 



where g~^ : R+ — > R+ is the inverse function of g that maps 
its range to its domain. Therefore, the total cost of multicast 
flow F is simply the sum of powers of all the hypearcs in the 
system. Note that the function g^^ is increasing and concave, 
and if h is convex then from Inequality (|2]l, g~^ oh increasing 
and convex. So for a given relay position Zr G C, the min-cost 
problem minimizing the total cost for setting up the multicast 
session (s, T) with a target flow F can be written as. 



Minimize 



(B) 
(6) 



subject to: F < F, < ^Fp , Vi G [l,ri], 

e:(P,Z?) 3 Fi. > 0,Vj G [l,T,]yi G (7) 

Constraint (|6]l makes sure that any destination ti £ T receives 
a minimum of flow F. Like in Program (A), we denote with 
the set D) the hyperarc rate and power constraints. 

Finally, define the point p*, that will be crucial in developing 
algorithms in later sections, as 

Zp. ^ a.i-gmm{ma.x{i^*h{D^^s), fJ-* '[na.x{h{D^^t,)))), (8) 
Zp tteT 

where, /i* = g{^) and ly* = g{v). An easy way to understand 
p* is that if ^* = V* = 1 then p* is the circumcenter of two 
or more nodes in the set A/^\{r}. Note that the program in 
Equation dH) is a convex program. Also, denote the optimal 
value of the objective function in Equation ([SJ as Dp-. 

Hereafter, we represent with (s, T, Z, 7) and (s, T, Z, 7, F) 
the joint relay positioning and flow optimization problem 
instances that maximizes the multicast flow and minimizes 
the total cost for a the target flow F, and with z*^ and Zp^ 
denote the optimal relay positions, respectively. 



III. Multicast Flow Properties And Reduction 

In this section we develop fundamental multicast flow prop- 
erties that govern the multicast flow in the wireless network 
hypergraphs that we consider in this paper. First, we briefly 
note the main hurdles in jointly optimizing the problem. For 
a given problem instance different relay positions can result 
in different hypergraphs, which makes the use of standard 
graph-based flow optimization algorithms difficult. Moreover, 
the hyperarc rate function can be non-convex itself. 

We will show that the joint problems {s,T,Z,j) and 
{s,T, Z,^, F) can be reduced to solving a sequence of two 
decoupled problems. The reduced problems are decoupled 
in the sense that the first problem is purely a geometric 
optimization problem and involves no flow optimization and 
vice versa for the second problem. At the same time, they 
are not entirely decoupled because the two problems need 
to be solved in succession and cannot be solved separately. 
Now we present a series of results that are fundamental to the 
reducibility of the joint problem. 

Proposition 1: The optimal relay positions z*^ and Zp^ lie 
inside the convex hull C. 

Refer Appendix lAl for the proof. Proposition [T] tells us that 
only the points inside the polygon C need to be considered. 
This brings us to the following fundamental theorem. 

Theorem 1 (Flow Concentration): Given Zj. G C; 

(i) the maximized multicast flow F* concentrates over at 
most two paths from s to the destination set T. 

(ii) for any target flow F G [0, F*] the min-cost multicast 
flow concentrates over at most two paths from s to T. 

The proof is detailed in Appendix |B] Theorem [T] is central 
to the two questions we aim to answer and reduces the 
complexity of joint optimization greatly by considering only 
two paths instead of many. Essentially, Theorem [T] tells that 
for a given relay position Zr G C, the multicast flow F 
must go only over the paths that span all the destination 
set T, i.e. set L. Furthermore, among the paths in L, the 
maximized multicast flow F* goes over only two paths, 
namely the path li = {(s, Ti), (r, T2)} that has the highest 
min-cut among all the paths through the relay r, and path 
I2 = {{s,ti, ..,tn) = (s,r)}, which is the biggest hyperarc 
from s spanning all the destination set T, where r G Ti and 
Ti U T2 = T. The same holds for the min-cost case for a 
given relay position Zr G C. Consequently, it is also true for 
the optimal relay positions z*^ and Zp^. Hereafter, we only 
need to consider the flow over paths li and I2 (corresponding 
to the relay position in consideration). 

A. Max-flow Problem - {s,T, Z,j) 

Assuming that the transmitted signal propagates omnidirec- 
tionally, we can geometrically represent the hyperarcs of the 
path li = {{s, Ti), (r, T-z)} by circles C^^ and C^^ centered at 
s and r with radii tTs = Dst^. and tt^ — F>rt^, (where Dgt^ = 
maxt,gTi(£'siJ and Drt^ = maxt^.gTa (-Drt^ )), respectively. 
Similarly, the path I2 = {{s,T)} can be represented by the 
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Fig. 2. The solid piecewise lineal' segment in examples (a) and (b) mai'ks 
the set of points r for different values of tt^ 6 (0, D^fj)- Each point r 
corresponds to 2*^ for some 7 G (0, 00). The piecewise linear segment breaks 
beyond the dashed circle as z\ G C|,^ . (a): E.g. C!^ with < TTa < Dst^, 
= arg min max(£)2-^jj , D^-^t^). Same goes for the example in (b). 

the union region of the two circles. Then using Theorem [T] 
Program (A) can be re-written as. 



Maximize ( min 



h{TTs) ' h{Trr) J h{Dst„) 



(C) 



where, P^_^ , and are the powers for hyperarcs of the 
paths li = {Cp_^, Cp^} and I2 = {C^}, respectively. The radii 
variables tts and tt^ correspond to path li for the relay position 
Zr G C such that Zr G C^^ and Z G Cu- 

Although Program (C) is reduced, it is still a non-convex 
optimization problem. The objective function is non-convex 
and different positions of the relay Zr ^ C result in different 
end node sets Ti and T2 for the hyperarcs of path li. 

On the other hand, we know that the relay position is 
sensitive only to the flow over path li. In addition, as there 
always exist a relay position Zr <E C such that the min-cut of 
path li is higher than that of path I2, then this also holds true 
for z*^. Therefore, optimizing the relay position to maximize 
the flow over path li results in global optimal relay position 
solving the original problem {s,T, Z,^). This motivates the 
decoupling of computation of optimal relay position from the 
flow maximization over the path li. 

Proposition 2: For a given problem instance (s, T, Z, 7), if 
g{v)h{Dsp') = Dp*, then z*^ = Zp*. 

Refer Appendix |C] for the detailed proof. At point p* , in 
general the following holds ^'"'tl > ^^'^1 (from Equa- 



tion (|8]i), thus making it naturally a good candidate for 
Proposition |2] essentially proves that if the relay is positioned 



'7f 



at p* and we get 



, and if maximizing the 



flow over the path li results in no spare source power (i.e. 
g{iy)h{Dsp*) = Dp,), then z* = Zp. and F* = 
Furthermore, the joint problem in Program (C) can be reduced 
to solving in sequence the computation of the optimal relay 
position p* by solving Equation (O and then calculating the 



max-flow F*. But this is not true when 



> 



We 



circle with radius Dst 



U CJp denotes 



cover this case in the section of algorithms. 

Let us now see the problem in a different way. Consider the 
radius tTs G (0,I?st„) and construct the hyperarc . Denote 



by T' = {tj e T\Dst - > TTs}, the set of destination nodes that 
He outside the hyperarc circle C^^ . Then compute the point r 
such that 

Zf = arg min( max (Z^r't, )) , 

and position the relay at r (here r is the point in C*^ such 
that the maximum among the distances to the nodes in the 
set T' from r is minimized). If Dgp < tTs, then we contract 
the hyperarc to C|, else we simply re-denote it with C^. 
Finally, we can construct the hyperarc Cj*^^ ( note that Z S 
Cy = Cp U ). The set TZ' of points r computed in this 
way for different values of tTs G {0,Dst^) are the optimal 
relay positions z*^ solving {s,T,Z,j) for some 7 G (0,cxd). 
Figure |2a) captures this interesting insight of the relationship 
between the points r and z*^. Note that the set TZ of points r 
is a discontinuous piecewise linear segment. 

B. Min-cost Problem (s, T, Z, 7, F) And Duality 

The min-cost problem (s, T, Z, 7, i^) can be written as 

Minimize + F^, + P|) (D) 



subject to: F < min 



h{ns) ' h{TTr) J h{DstJ 



Pt < Pk < ^• 



(9) 
(10) 



In the non-convex Program (D), the path li = {C-fi ' i 
correspond to the relay position Zr G C which is implicitly 
represented in the distance variables tTs and tt^. From Theo- 
rem [T] we know that paths li and I2 carry all the min-cost 
target multicast flow F. In this sub-section we refer the path 
li as the cheapest path for a unit flow among all the paths 
through r in i for given position of relay. 

Now, we claim that Zp^ G TZ. This is true because given the 
hyperarc Cf.^ of path li with optimal radius tt*, the second 
hyperarc C^^ must be centered at the point that minimizes the 
maximum among the distances to all the destination nodes not 
spanned by the hyperarc C^^ from itself, as this minimizes 
the cost over the hyperarc C^^. Therefore, Zp^ (like z*^) 
always lie on on the curve TZ. This observation motivates an 
interesting fundamental relationship between Zp^ and z*^. 

Theorem 2 (Max-flow/Min-cost Duality): For F G [0, F*], 



where 7 G [min(7, 7), max(7, 7)] and z*j_ 

Theorem |2] establishes the underlying duality relation be- 
tween the max-flow problem (s,r, ^,7) and the min-cost 
problem {s,T, Z,j, F) and says that the point Zp^ (or z^^) 



lies on the segment z* 



U 



(z^^ — z*^, respectively) of the 



curve TZ. Implying that the optimal relay position Zp^ solving 
the problem {s,T, Z,^, F) is also the optimal relay position 
solving the problem (s, T, Z, 7) for some 7. The proof of 



'7t 



Theorem |2] is presented in Appendix |D] 

However, the max-flow is not always reducible to a se- 
quence of decoupled problems. This is mainly due to the fact 



that the path I2 can be cheaper than path li for a unit flow 
corresponding to the optimal position z 



I.e. 



This information is not easy to get a priori. In contrast, we 
can safely assume that 



g-\hin:)) + g-'ihiK)) < g-'{h{D,tJ), 



(12) 



as almost all wireless network models that comply with our 
model result in the hyperarc cost function g^^{h{Duv^,^J) 
being the increasing convex function of distance Duv^^ th^t 
satisfy Inequality (fTZt . If Inequality (fT2l i holds, then similar to 
the Max-flow problem the joint optimal relay positioning and 
min-cost flow optimization problem in Program (D) can be 
reduced to a sequence of decoupled problems of computing 
the optimal relay position and then optimizing the hyperarc 
powers to achieve the min-cost flow F in the network using 
the similar arguments as in previous subsection. For a spe- 
cial of the min-cost problem {s,T, Z,^, F), we present the 
Min-cost Algorithm that sequentially solves and outputs the 
optimal relay position and powers to achieve the target flow 
F G [0,F*] in Section lEH 

IV. ALGORITHMS 

In this section we present the general max-flow and the 
special case min-cost algorithms that solve the sequence of 
decoupled problems. 

A. Max-flow Algorithm 
Input: Problem instance (s, T, Z, 7). 
1: Compute p*, if g{v)h{Dsp-') = g{ii)h{Dp-t^), output 



St 



= Zp*, F* = g(y)h{Dsp-') and quit, else go to 2. 



2: Construct the set T' = {t'^ G T\Dst'. < Dp^f ] = 
{t'l, ..,t'j,} (ordered in increasing distance from s) and 
compute . If -Dsf , < Dgp^^^^, declare 



^7t ~ ^Pt\t' 



and F* = g{iy)h{Dsp'^^^J and quit, else go to Step 3. 
3: Compute the points zjf and Z2, and maximized multicast 
flow Fl and F2, respectively. Declare before quitting. 



zl if F^ > F2*, 
z| if F^ < F2*. 



(11) Output: z*^ and F 



Fig. 3. Max-flow Algorithm. 

The Max-flow Algorithm in Figure |3] is a simple and 
non-iterative 3 step algorithm that outputs the optimal relay 
position and the maximized multicast flow. The first step is 
essentially Proposition |2] in case it is not satisfied the second 
step filters the redundant nodes that are too close to the source 
and can be ignored. If the conditions of first or second step 
are not met, then the third step divides the computation of z*^ 
into two regions of C and computes the optimal relay position 
zl and Z2 for these two regions and outputs the better one. 
The proof of optimality is provided in Appendix |E] 



B. Min-cost Algorithm 

In this subsection, we assume that the InequaHty ( fT2] ) is 
satisfied and the target flow F G [0, F*] goes over the path li 
(corresponding to the optimal relay position z*p^ only. Min- 
cost Algorithm in Figure IH unlike the Max-flow algorithm, is 
an iterative algorithm. In the first step the geometric feasibility 
region is constructed and in the second step this region is 
divided into at most n — 1 sub-regions. The optimal relay 
position is computed for all the sub-regions and the one 
minimizing the cost among them is declared global optimal. 
Computing the optimal relay position for the sub-regions is a 
simple geometric convex program that can be solved efficiently 
and the number of such iterations is upper bounded by n — 1. 
The proof of optimality is presented in Appendix |F] 

Input: Problem instance {s,T, Z,j, F) and C^. 
1: Compute p = argmin(/i(Dsp) + max {h{Dpt^))), and 

build the set f = {? e T\D^j; < Dsp}. If f ^ {0}, 

then recompute p = argmiii(/i(D5p)+ max {h{Dpt))), 

pec'r, teT\{f} 
calculate "i/p ~ h{Dsp) + Dpt^ and to go to Step 2. 
2: Build the set T = {i e T\{T,tn}\Dst > i^l.Dst < 
TT^} = {ti, ..,ti} (ordered in increasing distance from s), 
compute the points 

Pj = aTgmin{max{h{Dsp),h{D^j^_J) + max{h{Dpt))), 
pec] ' t^Tj 

and calculate the cost of unit flow "ij ~ h{Dsp) + 
h{masL{Dpt)) over the path I2 corresponding to the relay 

position pj, Vj e [1,/]. Declare 



if % < *m, 
if ^-ff > 



where = min (^-j), P| * = g-^{h{nl)F) and 
ie[i.n] ^ 

P£* = g-i(/i(7r;)F) and quit. 
Output: z*p^, P^^* and P^/. 



Fig. 4. Min-Cost Algorithm. 



V. Example: Low-SNR Achievable Network Model 

In this section we present an example from the interference 
delimited network model that was originally presented in [I]. 

A. Low-SNR Broadcast and MAC Channel Model 

Consider the AWGN Low-SNR (wideband) Broadcast 
Channel with a single source s and multiple destinations T = 
{ti,..,tn} (arranged in the order of increasing distance from 
s). From ||5] and IS), we know that the superposition coding 
is equivalent to time sharing, which is optimal. Implying that 
the broadcast communication from a single source to multiple 
receivers can be decomposed into communication over n hy- 
perarcs sharing the common source power Therefore, we get 
the set of hyperarcs Abe = {{s, ti), (s, ^1^2), ■•, {s, tit2-.tn)}- 



Similarly, in the Low-SNR (wideband) regime, interference 
becomes negligible with respect to noise, and all sources can 
achieve their point-to-point capacities analogous to Frequency 
Division Multiple Access (FDMA). In general, the MAC 
Channel consisting from n sources si,...,s„ transmitting to 
a common destination t can be interpreted as n point-to- 
point arcs each having point-to-point capacities. Thus, we 
get A,nac = .., (s„,<)}. Each hyperarc {s,ti..tj) e 

Abe U Amae IS associatcd with the rate function 



where a > 2 is the path loss exponent 



yje[i,n], 



(13) 



B. Low-SNR Achievable Hypergraph Model 

By concatenating the Low-SNR Broadcast Channel and 
MAC Channel models we obtain an Achievable Hypergraph 
Broadcast Model. For example the Broadcast Relay Channel 
consisting of a single source, n destinations and a relay. 
Although, the time sharing and FDMA are capacity achieving 
optimal schemes in the respective models, the Achievable 
Hypergraph Model is not necessarily capacity achieving. In 
contrast and more importantly for practical use, this model is 
easy to scale to larger and more complex networks. 

The above Low-SNR Achievable Hypergraph Model also 
incorporates fading [l]. The rate function in Equation (T3[ is 
linear in transmitter power and convex in hyperarc distance, 
hence the results from this paper can be directly applied. 

VI. CONCLUSION 

We present simple and efficient geometry based algorithms 
for solving joint relay positioning and flow (max-flow/min- 
cost) optimization problems for a fairly general class of 
hypergraphs. Any application that satisfies the hypergraph 
construction rules and can be modeled under the classical 
multicommodity framework can use the results presented here. 

As a part of future work it would be of interest to extend 
the work presented here to the general multicommodity setting 
where multiple sessions use a common relay. 
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Appendix A 
Proof of Proposition[T] 

Proof: Let the set of nodes M = {s, r, ti, .., t„} be placed 
on the 2-D Euclidean plane and C denote their convex hull 
polygon. Let us assume that the relay node r is placed outside 
the polygon C, i.e. <^C and c be the nearest point to r in the 
polygon C. Let the line segment joining and c be denoted 
as Zy — c. 

The rate over all the hyperarcs that either emanate from r 
or r is the farthest end node of the hyperarc, is relay position 
dependent. As the hyperarc rate is a decreasing function 
of distance, moving the relay closer to c on the segment 

— c decreases the distance between r and every point in 
the polygon C and thus to every node in the system. This 
implies that for a given power allocation for the relay position 
dependent hyperarcs the rate can be increased as the relay gets 
closer to the point c. Consequently, we can conclude that the 
optimal relay position z*^ maximizing the multicast flow F 
for the session (s, T) will lie in the convex hull polygon C. 

Similarly, all the relay position dependent hyperarcs will 
need lesser power to carry a given flow of value F as the 
relay r moves closer on the Une segment — c to the the 
point c. Implying, that for any target flow F the optimal relay 
position z*p^ will lie in C. This concludes the proof. ■ 

Appendix B 
Proof of Theorem[T] 

Before we formally prove Theorem [T] we need to establish 
some basic tools from convex analysis. 

Let / : R+ — > R+ be an increasing and convex function 
that maps a non-negative real input to a non-negative real 
output. Denote with /(x) the sub-derivative of j(x) at the 
point a; e R and let df{x) denote the complete set of sub- 
derivatives at point X. If the set df{x) is a singleton set, then 
f{x) = which is simply the derivative of the / at x; else 
there exist a finite interval df{x) between the left and right 
limits of / at x. In addition, let us also assume that /(O) < 0. 
Then the following proposition is true. 

Proposition 3: If / is any increasing convex function such 
that /(O) < then 



{xi 



^) e R" 



(14) 



Proof: As / is increasing over the real line, for xi < X2 
we have f{xi) < f{x2)- Also, as / is convex /(xi) < f{x2). 

Let the slopes of line joining the points (0, /(O)) and 
(xiJixi)), {xij{xi)) and {x2j{x2)) be given by, 

Sl = , S2 = , (15) 

X\ X2 - Xi 

where < xi < X2 are points on real line. From the 
Generalized Mean Value Theorem we know that there always 
exist a point c and c' between [0, xi] and [xi,a;2], such that 
/(c) ~ Sl and /(c') — S2, respectively. This, along with the 
fact that / is increasing and convex implies, 

/(O) < /(C) = Sl < /(Xl), /(Xl) < /(C') =S2< f{c'). 











» 








■3 









Fig. 5. (a): Increasing convex function / (possibly non-differentiable, e.g. at 
x'). (by. Orange, blue and red lines joining the points {/(O), 0) — (/(xi ), zi), 
(/(a;i),a;i) - {f{x2),X2) and {f{x3),xz) - (f(x2),X2), respectively. 

In general, given n points xi < ... < Xn with Si as the slope 
of line joining the points {xi, f{x,))-{xi+i, f{xi+i)) we get. 



Sl < S2 < ■■■ < s„ 



(16) 



Consider now the four points (0,/(0)), (xi,/(xi)), 
{x2,f{x2)) and (a;i2/(a;i2)), where xi2 = xi + X2- From 
Inequality [16] we get, 

/(Xi) - /(O) ^ f{x2) - fixi) ^ f{x,2) - f{x2) 



< 

Xi X2 - Xi 

Inequality [17] implies. 



Xi 



(17) 



(18) 



fi^l) + f{x2) - /(O) < f{xi2) - f{xi + X2). 

If /(O) = 0, then f{xi) + f{x2) < f{xi+X2). Using this fact 
it is straightforward to show that this also holds for /(O) < 0, 
for any {xi,X2) £ R'^^. 

Without loss of generality, assume that xi2 < X3, repeating 
the previous step of InequaUties [TtKTsI we get, 

/(X12) + /(.T3) < /(X123) ^ 

fixi) + f{x2) + /(.T3) < /(Xi +X2+ X3), 

where 0:123 = Xi+X2+x^. Similarly, repeating these n times 
we have 

n / n \ 

Y.f{x^)<fiY,^^]^ (X1,..,X„) eR"+, (19) 

1=1 \i=l / 

if ./(O) < 0. This proves the proposition. ■ 
Now let / and g be increasing convex functions satisfying 
Proposition [3] and define fi{xi) = \if{xi) and gi{yi) = 
X'ig{yi) as the 2n linear compositions of the function / and 
g fox i £ [l,n], where A; > 0, A- > 0, e and 
/(O) = g{Q) ~ 0. Then consider the following program. 

Maximize [ Fn = ^n\in{f,{xi),gi{yi)) \ (TIA) 



subject to: ^ Xi < /i, ^ yi < v, 



i=l 



x^ > 0, y^> 0, e [l,n] 



Denote the set 5** as the set of optimizers of Program (TIA). 
In addition, assume that 

min(Afc, Aj.) = max (min(Aj,A-)) , fcG[l,n]. (20) 

ie[l,n] 



Let us denote a set of points U* = {{xk,yk) = 
(0,..,0,M,0,..,0),(0,..,0,i^,0,..,0)} where, Xfc and are 
the vectors with /i and ly at the fc*'' place and all other elements 
0, respectively. Then the following proposition is true. 
Proposition 4: U* C S*. 

Proof: Consider program (TIA) and without loss of 
generality assume that 

min(Ai, A'l) < min(A2, A'2) < .... < miii(A„, A^), (21) 

and 

min(Afe, A^.) = min(AA;+i, A'j,+i) = .. = miii(A„, A^), (22) 

where Equation ( |22] | says that the last n~ k terms of Inequal- 
ity in} are equal. Let 

Pi — \ 

[gtivt) ^ Kgivi) if A, > A',. 

Then, from Proposition |3] for any < (e, e) G dom(/, g) we 
get. 



Proposition 5: (x^/,yk') C S'* . 

Proof: Let us assume that the inverse functions 

d g^"^ exists such that f^'^{\if{xi)) = Xi in addition to 

-1/ 



1=1 



which implies 



1=1 



min(/,(e),g,(e)) < min(/, (ne), ^^(ne)), 



th 



term 



as Pi is the limiting sub-term of the i 
mm{ f i{xi), gi{yi)) of function J'n for any i G [1,"-] 
and j G [k,n]. From Inequality ( l23T i. we can infer 
that simply maximizing nim{fj{xj), gj{yj)) (for any 
j G [k,n]) alone maximizes in Program (A), with all 
other terms attaining the value (except for j*^ term) 
because X,f{0) = X'J{0) = 0,Vi G [l,n\. Therefore, 
J-"* = min(/j(/i),5j(z/)) for any j G [k,n] is the maximum 
value of the function J'n in Program (TIA). 

Hence U* C 5*, where U* = {(xj,yj) = 
(0, .., 0, ^^, 0, .., 0), (0, .., 0, 1^, 0, .., 0)1.7 G [k, n]}. ■ 

Under the same setting as for Program (TIA), consider the 
following program. 



Minimize 



(TIB) 



V «=i / 

n 

subject to: C < ^'ncii't^iMx,), giivi)), 



i=l i=l 

x^>Q, >0, V^G[l,n], 



where C > is a given positive real number and such that 

C<min(/,(/.),g,;(z.)), V^G[l,n]. (24) 

Denote with S'*, the set of optimizers of Program (TIB). Then 
the following Proposition holds true. 



9i i-^igiUi)) — 2/i' for all i G [1, n\. Also assume that for any 
e > we get. 



f~MX,,f{xl,))+g-\Xk'9{yl,)) = xl, + 



Vk' 
mill (a;^ 

'ie[l,n] 



(25) 



where Xif{x\) = X[g{xl) = e,\/i G Since, / and g 

are increasing convex functions for non-negative input, their 
inverses /^^ and g^^ are increasing concave functions for 
non-negative input. From Proposition |3] we can deduce the 
reverse inequality for concave functions. 



(26) 



,4=1 



where /' : R+ 
From Inequality 

f~\Xk'finxl,)) 



< 



R+ is an increasing concave function, 
we get for any e > 

f .9-i(Afc-.g(ny|,)) 

(27) 



J2{f-\XJ{xl))+g-\Kg{yl))) 



(23) Lastly, from Inequality ( l24l i and (|27] | we conclude that 



K,,yl) C 5'*, where x^, = (0, .., 0, C, 0, 0) and y^, 

(o,..,o,c,o,...,o). ■ 

Note that we did not assume the differentiability of the 
functions / and g. 

The set of paths L, that span all the destinations will be 
central to the proof of Theorem [U so let us clarify some 
more notations. The path set L =- {/i,...,?,-} contains only 
one path from s to T that does not go through r, namely 
{(s, ti, .., r, .., t„) = (s,T)} and without loss of generality 
let us assume l-,- ~ is,T). All other paths go through r and 
consist of two hyperarcs, i.e. Ij ~ {{s,T(), {r,T2)} where 
r G Tf and U = T, Vj G [1, r - 1]. Then the path flow 
is given by 



'min(i?i,,i?: 



Rn 



where A] = ^ 



= min(A].g(P^,),g(A2p-p) 

ifjG[l,T-l], 

if j = T, 

V? G [1,T-1] and A, = 

^J^^ y where t{ £ T and t2 & T are the farthest nodes 

from s and r spanned by the hypearcs (s,Tf) and (r, Tj), 
respectively. 

Proof of Theorem 1: Consider the hypergraph G{M , A) 
for the given position of relay G C and the path based 
formulation of multicast max-flow and min-cost problems in 
Program (A) and (B), respectively. 

Since the hyperarcs are constructed in the order of increas- 
ing distance from the transmitter, there exist no two paths 
from the s to any destination ti G T that are edge disjoint. 
This implies that only the paths spanning all the destinations in 



the set L should to be considered, as sending any information 
over the paths that span a subset T' C T has to be resent over 
at least one path in the set L that spans aU the set T. This 
fact reduces Programs (A) and (B) to 



Max-flow 


Min-cost 


Maximize F = Fi (TIC) 

LeL 

subject to: < ^ F; , Vi 6 [1, n], 
Fi > 0, V/ e L. 


Minimize P (TID) 

subject to: F; < ^ F; , Vi G [1, n], 

F < Fi, Vi e [l.n], 
Fi > 0, V/ e L, 



where P = + P'h ) + -Rf i" Program (TIC), i.e. 

the sum of powers of all the nyperarcs of all the paths in L. 
Therefore, 



min Fi 

ie[l,ri] 



leL 



Without loss of generality let us assume that among all the 
paths from s through r to T the path Ik E L has the highest 
min-cut, i.e. min(A;^^ , A^^ ) > max_,g[i^_i] Xj,,Xf,. Then we 
get two scenarios, either 

min(A},A?) < .. < min(Ai,A^) < A^, (28) 

where, the last inequality of Inequalities ( |28] | says that the path 
(s, T) has the highest min-cut among all the paths in L. Then 
from Proposition |4] the multicast flow can be maximized by 
simply maximizing the flow over the path 1^ = (s, T), and 
since maximizing the flow over this path consumes all the 
source power fi the optimal multicast flow F* is given by 



F* = Fil = i?^* = XrOifi). 



(29) 



Otherwise if, 

min(A}, A^) < .. < A^ < .. < min(A^, A^), (30) 

then again by Proposition |4] maximizing the flow only over 
path Ik maximizes the multicast flow F in Program (TIC). 
Thus, we get 

Fil = min(i?-.,i?^t) = mm{Xlg{fi),Xlg{,y)). (31) 

Furthermore, if X\g{fi) < X\g{v), i.e. if the source has 
relatively more power than relay r, then the rest of the flow 
must be send over the path (s, T) as any other path through 
the relay (i.e. Ij where j ^ k and j G — 1]) cannot be 
used due to no spare power left with relay. This implies 



F* :=min(i?^1,,i?^P + i?f, 

= min(A^.g(Ai'), Xlg{iy)) + Xrg{n - n') 



(32) 



where Xlg{n') = Xlg{iy) and i?f, i?f,(P| ^ l-i - n') = 
XrgilJ. — m')- Thus, all the maximized multicast flow F* goes 
over at most two paths, l^. and l^-. Integrating Equations d3TT i 
and ( [32] i. we get 



if A^ = max (min(Aj, A^)), 

i6[i.T] 

F* = <! min(i?^^ , + i?l 

if min(A^,A^) = max^(min(A], A^)). 



(33) 



From Equation ( |33] |. we conclude that for any given position 
of relay Zr G C, the optimal multicast max-flow F* goes over 
at most two paths namely Ik and Ir- Consequently, this also 
holds true at z*^. 

For the case of multicast min-cost Program (TID) for the 
target flow F £ [0, F*], without loss of generality let us 
assume that 



Afe'+A^ = max (Aj+A^). 

je[i,r-i] 



(34) 



From Equation (|5]i cost of sending the flow e > over the 
path Ik' is given by 



^'.^-^ (aL 



where g ^ is the inverse function of the power function g and 
is increasing and concave. From Equation ( |34] | we get 



min {P', +P; 



(36) 



For a given position of relay G C, then clearly 
min(A^/, A^., ) = min(Aj.,A^), i.e. the cheapest path through 
the relay is that path with the highest min-cut. This is true 
because 



min(A5',A2)= max (miii(Aj, A^)) 

je[i,T-i] 



(37) 



= mm 



A^ X\ je[i7r-i] \ A^ Xi 



(2) P^., 

(3) P^,, 

(4) pL 



if A| = maxjg[i^_i] (A2) and this can be safely assumed for 
the path with the highest min-cut. 

From Equation ( l36t and Proposition|5l we infer that ZJ. is the 
cheapest path for a unit flow among all the paths 1^ G [1 , r — 1] . 
Moreover from Equation ( |37] |. paths Ik' (= Ik) and l^- can carry 
any target multicast flow F G [0, F*\. So we get four cases 

(1) P^„^^ + P;,,^^ < Pf , and miii(Ai,, A2,) > A„ 
P^y ^ < PI, and min(Ai,, A2,) < A,, 
P;V^>P|,, and min(Ai,,A2,)>A., 
PL,/ > Pf , and min(Ai,, A^,) > A^. 

In all of the above cases, the target flow F G [0, F*] 
flows over the paths Ik' and Ir only. Thus, we conclude that 
for any relay position z,. G C the optimal min-cost target 
multicast flow F flows over at most two paths Ik' and Ir, and 
consequently also at the relay position z|,^. This completes 
the proof. ■ 

Appendix C 
Proof of Proposition!!] 

Proof: From Theorem [T] we know that the maximized 
multicast flow goes over at most two paths, namely path li 
having the highest min-cut among the paths through r and 
path I2 that spans all the nodes in the system. Moreover, there 
always exist at least one relay position such that the min-cut of 
the path li is at least as that of path I2, implying that this also 



holds at the optimal relay position z*^ solving {s,T, Z,j). 
This is also true for point p* because at p* 

ma.x{h{Dsp*),h{Dp.tJ) < h{DstJ, 

that comes from its definition in Equation 

Positioning the relay at p* will render the highest min-cut 
of path li compared to that for any other position. This is true 
from the definition of point p* itself. If at point p* we have 



g{iy)h{D, 



Dp, = min,g[i „] g{^)h{Dp*ti), then 



F*, = min 



where, Trf = Dgp' and tt^ 



HG l,n 



From our assumption at point p* the maximized flow F* 
consumes all the source and the relay powers. Since we only 
consider those positions of relay at which the min-cut of path 
li is higher compared to path I2, positioning the relay at any 
point p E C such that tt^ > ttJ only renders decreased 
maximum rate over the hyperarc C^_^ of the path li. Implying 
that F* < F*, , even though there might be some relay power 
left. 

On the other hand, positioning the relay at point p such that 
Trf < 7rf , increases the maximum rate over the hyperarc C|.^ , 
as 



gjfj-) 

h{Dsp[ 



> 



.9(m) 
h{D,p,)- 



Moreover, we get ttJ? > tt^ , as moving away in any 
direction from point p* increases max^gT-^ (Dp.^^ ). Therefore 
the multicast flow is at this point is given by. 



p* 



r ^ a(u-a-^ ( sMM-Il]] 



h{n?) h{DstJ 

where the first term on the right hand side is the flow over the 
path 1 1 = {C^_^ , CJp^} that is limited the hyperarc C-J^^ and the 
second term is the flow over the path I2 = {C^} such that 
Fp is achived by maximizing the flow over the paths li and 
I2 successively. 

As a Corollary of Proposition |3] it can be seen that 

aigixi) + a2g{x2) < a^g{xi + X2), 

where g is an increasing convex function, > for i G [1, 3] 
are some constants such that 03 > max(ai,a2). This implies 
that for any source power e > 0, the flow over the path li 
corresponding to the relay position p* will always be larger 
than the sum flow over the paths li and I2 corresponding to 
the relay position p. Therefore, for any such relay position p, 
Fp < Fp,. This proves the proposition. ■ 

Appendix D 
Proof of Theorem[2] 

Refer Figure |2] as a reference example along with the proof. 
Proof: Let the optimal relay positions z\ and Zp. be 



respectively. Then the hypearcs of the path f| = {Cf,^, CX,^} 

-J 

and = {C^^,CZ,j^} can be constructed simply by forming 
the first hyperarcs Cf,^ and C^^ with radii ttJ* = Dsz' and 
7rf* = Dszi.^, respectively. Here, the paths l\ and li represent 
the path li corresponding to the relay positions z*^ and 
respectively. Compute the points 

z^t — argmax (max(Z?rt)), 
rGC%,zeCu teft 

n 

Zj.1 = argmax (max(Z?rt)), 
rGC= I ,zeCu tef J- 



where f ^ = {< e T\Dst, > 4* = -D..;,} and f^ = {te 

If the points and are not the same as z*^ and Zp^, 
respectively, then this contradicts the optimality of the two 
points z*^ and Zp^. This is true because the only way to either 
maximize the rate or minimize the cost over the hyperarcs C!^^ 

and CJ^^ is to compute the points inside the hyperarcs Cf,^ 
^2 ^ 1 

and Cf,4 that minimize the maximum among the distances to 

all the destination nodes outside these hyperarcs from itself, 
respectively. Therefore, the optimal relay positions z*^ and 
Zp^ solving the problems (s, T, Z, 7) and (s, T, Z, 7, F), are 



the points of type T on the curve 7^. Hence, Zp 



for 



some 7' G (0, 00). 

Now let us consider that the position z^^ that minimizes 
the cost of unit flow (normalized, if the flow values are less 
than unity) from s to T, and z*^ ~ z^^ for some 7 € (0, 00). 
Without loss of generality, let us assume that z^^ is situated 
on the right of z*^ on the segment TZ. This implies that the 
rate over the source hyperarc Rj,^"''^ is the limiting term in 

^::.^=min(i?;;;^^i?;f^), 

where FZ* is the maximized flow over the path /T"'^ cor- 

'1 

responding to the relay position z^^. Furthermore, the only 
way to increase the min-cut of the path li is to position the 
relay on the left of (closer to s and z*^) on the segment 
TZ, as positioning the relay further on the right of z^^ on the 
segment TZ will not only increase the cost of unit flow but 
will also decrease the min-cut of the path li. Therefore, for 
any F e [0, F*], z*p^ e z*^ - zl^, where z^ - zl^ is the sub- 
segment of the piecewise linear segment TZ joining the point 
z*^ and z^^. The same argument holds for the case when the 
point z*^ lies on the right of zi-^. Thus, we can conclude that 



for some 7 G [min(7, 7), max(7, 7)], and VF G [0,F*]. 

Now suppose that for some 7 G 7 — 7 such that ^ — jjj 
and [fj! ,v') < (Mj '^)' the optimal relay position maximizing 
the multicast flow F is given by z~^ G z*^ — z^^ and the 
given that solve the problems (s, T, Z, 7) and (s, T, Z, 7>), maximized multicast flow is given by F*. Implying that, we 
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Fig. 6. Consider i?^^^^ = t^^jj^;^ = 
in the examples, (a): \T\ = 3 example with 



with g(^t) = 



= 1 
and 



the path li = {C^^, CJp^} (corresponding to p*) carries all the maximized 
flow F*. Therefore, z*^ = Zp* . (b): |T| = 4 example showing the path 

= {Cr^'C'TaJ' corresponding to p* but ^^^^^^^ > ■(D~r71^- 
nodes in T' = {ti , t2} can be ignored, as the hyperarc CJ^ of the path Zi = 
{C|,^ , CJp^ } corresponding to p'^^j,, spans the set T'. Thus z*^ = Zp» ^ . 

need at least the total source and relay power of value fi' + v' 
to achieve the multicast flow F*. . Thus, we can conclude that 
z-^ = z*p* ^. Hence, for any 767 — 7, there exist a flow 

value F^*. S [0, F*\ such that = z*p,^ ^. This completes 

the proof of Theorem |2l ■ 

Appendix E 

Proof of optimality of Max-flow Algorithm 

We divide the space of max-flow problem in three cate- 
gories, each for a step in the Max-flow Algorithm in Fig- 
ure [3] Proving the optimality for each category will prove the 
optimality of the Max-flow Algorithm as any given problem 
instance (s,T, Z,7) will fall in one of the three categories. 
The first two categories (corresponding to the first two steps 
of the Max-flow Algorithm) deal with the instances when at 
the optimal relay position z*_j. all the maximized multicast flow 
F* goes over the single path l\. 

Proof of Step 1: Refer Appendix |C] for the proof of 
Proposition 12] and Figure |6j a) for an example. ■ 

Proof of Step 2: Now assume that g{v)h{Dsp*) < 
g{fi)h{Dp*t^). This implies that there exist at least one 
destination node that is closer to s than to p*, thus we can build 
the ordered set T' = {t'^ G T\Dst'. < Dp.tr)} = {t[, ..,t'^,} 
in increasing distance from s. 

We can then recompute the point p* for the destination 
nodes in the set T\T' and denote it by p^^-j,,. If we get, 
g{v)h{Dspr^^^,) = g{^i)h{Dp*^^^^tJ, using Proposition E] 
we infer that the point p^yj., is the optimal relay position 
maximizing the mutlicast flow from s to the set of destinations 
in the set T\T' . Furthermore, if Dgt',, < ^sp* then all the 
nodes in the set T\T' are spanned by the hyperarc Cf,^ with 
radii tt^ = Dsp*^^^, corresponding to the relay position p'^yj,,. 
Implying that the z*^ = Zp-^^^^ . ■ 

Step 2 essentially gets rid of redundant destinations for 
computing the point p* that are close enough to the source 
s. Refer Figure |6] for an example. 

In contrast, if at point p1^\^rp, we get g{v)h{Dsp*^^^^) < 
g{fi)h{Dp^^^^t„), then we can divide the problem of optimal 



relay position within two regions in C. Since the relay position 
maximizing the multicast flow in the two regions is going to be 
unique, we can compare the two results and declare the global 
optimal position z*^ solving the problem instance (s, T, Z, 7). 

First region is the interior of the circle intCp centered at s 

j' 

with radius Dgf and the second region as the area C\intC^, 
J j' 
which is the rest of region in C minus the first region. Let us 

denote the optimal relay positions inside the first and second 

regions with z* and Zj, respectively. At first, we present the 

following two propositions that will come in handy for the 

proof of optimality of Step 3. Refer Figure |7] with g{P) = P 

and h{D) = for an example. 

Proposition 6: For any |r| = 2 case max-flow problem in- 
stance such that g{v)h{Dsp*) < maxjg[i 2](.9(M)^(^p*ti)) = 
Dp-, zl 6 intC('*^ lies on the line segment p* and pi2. 

Proof: Consider the two destination node case such that 
g{v)h{Dsp') < max,^[l^2]i9{^J')HDp^u)), where 

p* ^ a.rgmm{ma.x{g{i^)h{Dsp), g{fi)h{Dpt^), g{n)h{Dpt2)))- 
This implies that Dgt', > D 



sp 



T\T' 



where 



p*, = aTgmm{ma.x{g{iy)h{Dsp),g{fJ.)h{Dpt^))), 

p£C 

and T' — {ti}. Moreover, as a consequence of Theorem [l] 
only those points need to be considered such that positioning 
the relay at this point gives higher min-cut for the path li than 
that of path I2, inside the region C D intCj^ . 

The circle C^, is the circle Cf^ centered at s with radius 

Dsti - Let the perpendicular bisector of the nodes ti and t2 be 
denoted as ±12, and its intersection point with the segment 
s — t2 by pi2. Now consider any point p e C fl intCj''^ in the 
halfplane containing the node t2 of the bisector ±12 and let 
p' denote the closest point on ±12 to p. Then clearly, Dsp' < 



D, 



h{Dsp') < h{Dsp) and Dp,t, = Dp,t,_ < I?; 



h{Dp't.-^) < h{Dpt^). This implies that the min-cut of the path 
li reduces when the relay is positioned at the point p compared 
to that at point p' . Therefore, for any relay position in the 
halfplane (of _Li2) containing t2, there always exist a position 
on the segment p* — pi2 that is a better candidate for the 
optimal relay position maximizing the multicast flow. 

Now consider a point p G C D intCf^ in the halfplane (of 
-L12) containing the node ti such that Dsp < Dsp', then we get 
Dpt2 > Dp>t2- Note that, positions p for which Dsp > Dsp' 
are uninteresting due to the fact that not only the min-cut of the 
path li decreases but also that the maximum possible rate over 
the hyperarc Cf.^ decreases compared to the relay position p' 
(corresponding to p). Then the maximized multicast flow for 
the relay position p and p' such that Dsp 



< Dsp' is given by 



" KDpt,) 



F*, = 
P h{D. 



h{D„,) 



h{Dst,) 



^12 




Fig. 7. |T| = 2 node example with T' = {ti}. Note that, Dsp < D^p, and 
Dpt2 < Dpif^. Therefore, the rate over the paths = {(s,p), (p,tit2)} 

and = {{s, p'), {p' ,tit2)} is limited by the hyperarcs {p,tit2) and 
(p',tit2), respectively. 

respectively. The first term in the above equations is the flow 
= m.m{R^^ , R^^) = R^^ over the path li which is limited 
by the rate RTp^ over the hyperarc C^^ and the second term 
is the flow over the path I2 which is the function of total 
source power /i minus the power used over the hyperarc Cf.^ . 
As Dp't2 < Dpt^, we have ) < ^(0"'! ) ■ Furthermore, 

from triangle inequality we get Dgp' — Dgp < Dpn^ — Dpt^, 
implying that 

g(£) _ g(e) < .g(£) _ g(e) 
h{Dsp) h{Dsp') - hiDp^t,) h[Dpt,y 

for any P^^ = = e > 0. Therefore, we get more flow over 
the path li corresponding to the relay position p' compared to 
that of relay position p. Lastly, since the min-cut of the path 
I2 is strictly smaller than that of path li, it can be concluded 
that Fp < F*, . Hence, for any point p e C Ci intQ^ there 
exist a point p' on the segment p* — pi2 of the perpendicular 
bisector ±12, such that F* < F*,, hence zl e p* - pi2- ■ 
Note that if h is an strictly increasing function of distance, 
then the following inequality would hold strictly F* < F*,. 
For |r[ > 2, Proposition |6] can be generalized in the following 
way. For simplicity, let us assume that the line £st„ passing 
through s and t„ is horizontal and point p* lies above £st„ 
(ref. Figure |7] for example). Now, let the set of perpendicular 
bisectors for each pair of nodes in T be denoted by ±= 
{_Li2, -L13, .., _L„_i„}, where _La6 denotes the perpendicular 
bisector of the nodes ta £T and tt G T, and | _L | = „|(,"i2)! ■ 
Most of the bisectors ±ab^^ will intersect with the line 
£stn ™d let Zab denote the angle between the point s, the 
point of intersection of Lab and £st^, and any point on the 
bisector _La6 above the hne £st„- The closest point on the 
perpendicular bisector Lab to s is denoted by p^^. In addition, 
let ±D ± = {_Lxrn -L2TH -LjTm} be the subset of bisectors 
_Lj„ of the nodes tj G T and t„ (the farthest node from s 



in the system) for j G [1,?ti], such that there exist a segment 
P]n ^ Q]n of -Ljn C SO that the farthest nodes from any 
point p G Pjn — q-j„ are the nodes tj and i„ themselves. For 
example, in Figure [T] the segment p* — pi2 of bisector ±12. 
Finally, without loss of generality we assume that ^[„G i be 
the perpendicular bisector such that the point pf^ hes on the 
segment - qj^. 

Proposition 7: For \T\ ~ n > 2 max-flow problem instance 
such that g{v)li{Dsp') < Dp. = ma.x,^[i^n]ig{ti')h{Dp*u)), 
z* G C n intCj**, ^ lies on the piecewise linear segment {p* — 

QUn Qln ~ '?2jii Ql-ln ^ tin) ■ 

Proof: Consider the region C H intC*/ , for determining 
the best possible relay position z* maximizing the multicast 
flow F. p* is already a good reference point in C fl intCj, ^ . 
From Proposition |6l we know that position p such that 
Dsp < Dsp', i.e. the rate over the hyperarc C^_^ increases 
for a given power is an interesting position in terms of 
finding the point z*, thus we will only consider such directions 
from p*. In other words, only those bisectors -Lj„G 1. need 
to be considered that intersect £st„ on the segment s — i„ and 
make an obtuse angle Zj„. 

Let tj be the limiting node in determining p*, i.e. 
g{li)h{Dp*ti) — Dp' whose bisector makes the largest obtuse 
angle Zj„ with 1st,-, (e-g- node ti in Figure [8ja)). Also, 
and tn are the farthest limiting nodes from s in determining 
the point p* . This implies, that there exist a segment p* — qi^^ 
on the bisector ±j„ (towards s) such that ti and t„ are the 
farthest nodes from any point p G p* — qi^- If pf^ G p* — qin, 
then using Proposition |6] for any position p G C fl C^, and the 

closest point pj^^ to p on -Li„, we get three cases 



either 




< ^-^PL 


and 


Dpt„ 


> ^P'lJ- 


or 




> Dsp'^^^ 


and 


Dpu 




or 


^ sp 


> ^■'^PL 


and 


Dpu 


< Dp>_j,_ 



Then if p^^ G Pi„ — qin, using Proposition |6] we can deduce 
that z*^ & P* — qin- All other points p G C n intC|, , such 
that pj^ lies outside the segment p* — gj„ can be ignored. 

On the other hand, suppose that pj^ ^ pj„ — qi„. Then, 
there exist another bisector J-2n (of the nodes t2 and t„) that 
intersects Zj„, say at point gj„ and contains a segment qi„ — 
q2n G CnintCf, , such that for any point p G qi„ —q2n, ^2 and 
tn are the farthest nodes in the system. Using proposition |6] 
again, we can infer that for all the points p G CCiintCp ^ such 
that Dpp>_ < DppL , positioning the relay at pj^ will render 
F* < F*, . Similarly, if Dpp'_ > Dpp'_ , then positioning the 
relay at p£ will render F* < F*,_^ . Finally, if p|^ G qin-qin, 
then we can conclude that z*^ g"(pi„ - qi„, qi„ - q2n)- 

Generalizing to case of / such bisectors such that pf G 

^ ^ In 

Qi-in - Qin^ we conclude that zl G (pi„ - qin, qin - 
l2m ■■Tli-in ~ lin)- The Same argument would suffice if p* 
lies below the hne £st„ and this completes the proof. ■ 

Now we are ready to prove the optimality of Step 3 of the 
Max-flow Algorithm. 
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Fig. 8. |T| = 3 example with T' = {ti}. (a): ±= {±12, ±13, -L23}. (b): 
_L = {-Lj3, -L23}, where tj = ti and = t2- ^^'^ is the angle bisector of 
Zp*pj^ps dividing the region CnC"_^ into the two halfplanes of points that 
are closer to segments PJ3 — gj3 and qj^ — 523- 

Proof of Step 3; Following Step 2, assume that not all 
the nodes in the set T' lie inside the hyperarc Cf.^ with radius 



D 



then by reforming the hyperarc C^_^ with 



radius tTs = Z^gj' we can compute the point 

p = argmin( max (Dpt^)). 

p€Ci, *iGT\{T'} 



(38) 



The point p will always lie on the circumference of the 
hyperarc circle Cf^^. Positioning the relay r at p gives the 
hyperarc C^^ with radius tt^ ~ Dpt^, thus rendering the path 
li = {Cf-^jCJ-^} for the relay position p. Since i^sp* ^, < 
Dgt' = Dsn, this implies g{i^)h{Dsp) > g{fi)h{Dpt^). Thus 

maximizing the rate over the path li will alone maximize the 
multicast flow F for the position p, giving 



= mm — -, — r 

P \h{Dsp)' h{DptJ 



(39) 



For any relay position p with Zp G C\mtC^_^ the min-cut of 
the path li reduces compared to the position p simply because 
Dsp < Dsp, implying that the maximized multicast rate F* < 
F?. Therefore, 4 = z|. 

For all the positions inside the circle C^, , from Proposi- 
tion IZl we know fliat e {p* -qi,„ qi„ - q^n ,..,qi^i„- qj^ ) , 
where {p* - qin,qin - 92™, 'j'r-in - Un) is piecewise 
linear segment made up of the sub-segments of perpendicular 
bisectors ± = {_Lj„, .., _L^„}, where from any point on the 
bisector ±fc„, t^. and i„ are the farthest nodes in the system, 
for all k G [1, Let us re-denote this piecewise linear segment 
by s = (sj„, .., s^„), where denotes the sub-segment 
Qk-in- Qknfor k e[l,n]. 

I : is easy to compute the equation of each bisector in ±. 
Let the equations be given by 



y = mix + ci, Vfc G [1, /], 



(40) 



where y = tojx + cj is the equation of the bisector and 
the segment qfc_i„ — q^.^ t)^ represented by limiting the 
values of a; G [xqi^_^, Xg^], where Zg-^_^^ = {xq-,_^^,yq^,_J 
are the coordinates of the point qfc_i„ in the plane. Moreover, 
as we know that for any point G qk-in ^ Ikn^ '^^e limiting 
hyperarc is C^^, implying that the maximized multicast flow 
Fp^ can be achieved simply by maximizing the flows over the 



paths li and I2 in succession, we get 



F* 

Pk 



HDp-,t-J 



HDstJ 



where the only variables are Dspj_ and Dp-^tj^ as pj. G qk-in ^ 
qj.^-^. The variables Dgpj, and Dpj^tj. further are functions of 
coordinates {xpj,,yp^) of point pj,. Using the Equation ( l40b . 
we can rewrite F*^ = F^{xj.) as a function of single variable 
Xf., where xj. G [xq-j^_^,Xq-^\. Then the optimal relay position 

maximizing the multicast flow F in the region C n mtC^, is 

3' 

given by 2* = {x\,yl), where 



x\ = argmaxf max Fr (xj.)), 
fcG[i,r] 



(41) 



where x^, G [x^^.^^ , Xg^J, Vfc G [1,^] and yl = mj.Xj^ + Cj.. 

Finally, comparing the values of F^ and F2*, the optimal 
relay position solving {s,T,Z,j) is given by 



"■it 



if F,* > Fi. 

if Fi* < f;. 



This completes the proof of optimality of Step 3 of Max-flow 
Algorithm. ■ 

Remark on solving Equation X4J\ : The single variable 
function F|(a;^) is non-convex and smooth over the domain 
Xf. G [xq^^^,Xqj^],\/k G [1,^]. Moreover it can be proven that 
there exist a single stationary point (maxima) of the function 
in the domain [.Tg-^,.Tq-]. Therefore, using gradient based 
approach the global maxima can be achieved, implying that 
non-convexity is not a hinderance. 

Appendix F 

Proof of optimality of Min-cost Algorithm 

Assume that for a given problem instance (s, T, Z, 7, F) the 
Inequality (fTSl i holds and at the optimal relay position Zp^ all 
the target flow F goes over path li only. Refer Figures l9lfT0l 
for example. 

Proof of Step 2: Assume that at the optimal relay 
position, all the min-cost multicast flow will go over path I2 
only. This implies that the intersection region = C"'' flC"*" 



-1 I g(tJ-) 

F 



and 



of the circles and C"" with radii tt^ = h 

TTj^ = (^F^V respectively, contains at least one point, 
and this assures feasibility. Positioning the relay at any point 
inside Cp will fetch the multicast flow of value F over the 
path l2- Implying that the relay position in minimizing the 
cost of unit flow also minimizes the cost of flow F and is thus 
the optimal relay position solving (s,T, Z,7,F). 
Now consider the point 



argmiii(/i(F)5p) 

peC'r, 



max{h{Dpu))). (42) 

i£[l,n] 



It can be seen that all the nodes in the set T = {t E T\D^ij~ < 
Dsp} can be ignored as they are close enough to the source 



Fig. 9. Consider the \T\ = 4, node system, witli = n C'*'i and 
tlie point p tliat is tlie optimal relay position minimizing the cost of unit flow 
over the path h in the region n C-. 

node to be spanned by the source hyperarc C^^ of the path 
li. Therefore, we can recompute 

p = arginin{h{Dsp) + max {h{Dpt^))), 

peC'r, teT\{T} 

and get rid of unnecessary bias. Denote the cost of unit flow 
= h{Dsp) + h{Dpt^ ) corresponding to the relay position 

P- 

Now consider the region C'f^r\C~ , where the circle C'^ has 
the radius Dsp- p is the optimal relay position in the region 
minimizing the cost of flow F, as all the destination nodes 
inside the circle are not the limiting nodes determining 
the position p. But the point p is not necessarily the global 
optimum for the whole region C'^. So we break the problem of 
finding the global optimal relay position minimizing the cost of 
multicast flow F, into finding the optimal relay position among 
disjoint regions of and then compare them to declare the 
global optimal relay position. 

Denote the set of no_des T = {t e T\{f,tn}\Dst > 
ttP, Dst < T^s} = {ii, --jti}, where 7rf = Dsp- Consider find- 
ing the optimal relay position pi in the region C'(^ri{C2\Ci) 
that minimizes the cost of unit flow over the path li, where 
circles C'l and C2 are centered at s with radii ZJ^j^ and Dgj^, 
respectively. Then the problem can be stated as. 

Pi = argmin(max(/i(Dsp), h(D^j )) + ma_x{/i{Dpt))), (43) 

where the set Ti = {t e T\Dst > D^jJ consists of 

destination nodes lying outside the circle Ci- The Program 
in Equation ( |43] | outputs for the optimal relay position mini- 
mizing the cost of unit flow over the hyperarc li for the relay 
position in the region fl {C2\Ci). Alhtough, the region 
Cp n {C2\Ci) is non-convex, the optimization program in 
Equation ( |43] ) is convex and easy to solve. Similarly, for all 
the nodes in the set tj e T, we can compute 

Pj = argmin(max(ft,{Dsp), /i(D^j . )) + ma_x(/i(Dpt))), (44) 

and the cost of unit flow at the optimal relay positions pj in 
the disjoint region Cp fl {Cj_f_i\Cj) can be calculated, thus 
denote the vector = {'I'l, .., ^f/}. Furthermore, generating 
the set needs < |T| = ?7 iterations. Finally, we get 



Fig. 10. |T| = 4 example of Figure|9]is considered with T = {ti,t2}. (a): 
pi on the circumference of dotted circle Ci is shown, which is the optimal 
relay position solving (s, T, Z, 7, _F) in the region C2\Ci. (b) For the region 
C"'\C2 and the optimal point p2 is shown on the circumference of C'2- 

izp if 
""^^'[zp^^ if%>l^^, 

where, = min^gfi^;] ('I'fe). 

Note that, the dividing of the region (7^ into the disjoint 
non-convex regions C'^\intCj lets capture the idea of po- 
sitioning the relay anywhere in the region (7p such that the 
source hyperarc C^_^ of the path li spans all the destination 
nodes that are closer to source than relay and this ensures the 
global optimality. This completes the Proof. ■ 



